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Abstract 

In this work we study the Dirac quasinormal modes of higher dimensional charged 
black holes. Higher dimensional Reissner-Nordstrom type black holes as well as charged 
black holes in Einstein Gauss-Bonnet theories are studied for Fermionic perturbations 
using WKB methods. A comparative study of the quasinormal modes in the two different 
theories of gravity has been performed. The beahviour of the frequencies with the variation 
of black hole parameters as well as with the variation of spacetime dimensions are done. 
We also study the large multipole number limit of the black hole potential in order to 
look for an analytic expression for the frequencies. 
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1 Introduction 



The dynamical evolution of perturbation in a black hole background [1^ 2} mainly consists of 
three stages - the first, an initial wave burst which depends on the original form of the field 
perturbation, second, the damped oscillation whose frequency and damping depends entirely 
on the background spacetime and not on the perturbing field. The final stage involves a power 
law tail behaviour of the waves at very late times. It is the second stage which is of most 
importance in studying black hole physics. This is characterized by a discrete set of frequencies 
termed quasinormal frequencies. On perturbing the black hole background, either through 
external fields or via metric perturbation, the black hole system responds by emitting damped 
oscillations [H [2]. The frequencies consists of a real part which represents the frequency of 
the oscillations and an imaginary part representing the damping. For a large class of black 
holes, the equations governing the perturbations can be cast into the Schrodinger like wave 
equation. For asymptotically flat space-times, the quasinormal modes (QNMs) are solutions 
of the corresponding wave equation with complex frequencies which are purely ingoing at the 
horizon and purely outgoing at spatial infinity [SI H]. One of the most important interest in 
studying QNMs arises from the fact that there is a possibility of detecting them in gravitational 
wave detectors |5l [6]. Apart from that, QNM's have been found to be a useful probe of the 
underlying space-time geometry. QN frequencies carry unique information about black hole 
parameters [7j. It has also been shown that quasinormal modes in Anti-de Sitter space-time 
appear naturally in the description of the corresponding dual conformal field theories living 
on the boundary [HI [9]. Therefore, many researchers found it interesting to study QNMs 
in the background of asymptotically AdS black holes [TOl [TT] . The QNMs of black hole in 
asymptotically de Sitter spaces were also studied accordingly [12], since there were evidences 
that our universe is governed by positive cosmological constant [13]. Let us also mention here 
for the sake of completeness that in spite of their classical origin, it has been suggested that 
QNM's might provide a glimpse into the quantum nature of the black hole [13]. 

Black hole QNMs have been studied since long time. A lot of studies were made for pertur- 
bations of black hole background with fields which are of integer spins (scalar, electromagnetic 
and gravitational). Compared to this, the study of Dirac QNMs in black hole backgrounds is 
rather very limited. The Dirac QNM for (2+1) dimensional BTZ black hole in AdS space was 
studied both numerically and analytically in [SJ [TT] . Studying Dirac QNMs in four dimensions 
started with the work of Cho [15] who first studied massive and massless Dirac QNMs in a 
Schwarzschild black hole background using the WKB approach [¥l] B2j followed by Jing [16] , 
who modified the approach using continued fraction [17] and Hill determinant method [18]. For 
the massless Dirac field, it was shown that the real part of the QN frequency increases with the 
multipole number I while the imaginary part increases with the overtone numbers n. For the 
massive Dirac field the real part of the QN frequency increases with the mass of the field, while 
the damping decreases, which implies that it would be possible to detect QN frequencies due 
to perturbations of black hole background with massive Dirac fields, since fields with higher 
masses will decay slowly. The high overtones of Dirac perturbations of a Schwarzschild black 
hole were studied in [T^]. Then there were lots of similar works on the Dirac perturbations of 
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four dimensional black holes both in asymptotically flat and non-flat black hole backgrounds 
[201 [21]. 

All the abovementioned studies of Dirac QNMs were done in the four dimensional black 
hole backgrounds while, the study of Dirac perturbations in higher dimensional black hole 
backgrounds are even more limited. String theory requires the existence of extra spatial dimen- 
sions and for a long time it was thought that the only possible way to think of these as extra 
spatial dimensions tightly curled with the radius of curvature around string scale. However 
later it was realized that these extra spatial dimensions need not be of the order of string scale 
but they can be as large as a few millimeters [22l [23l |2lj. An important feature of the large 
extra dimensions is that it implies that the fundamental Planck scale is much lower than the 
four dimensional Planck scale, in fact it might be of the order of a few TeV. This provides a 
tool for looking into stringy effects which might become observable at the LHC One of the 
most interesting effects will be the production of microscopic black holes at the LHC. Since 
gravity will be sensitive to macroscopic extra dimensions, these black holes produced at the 
colliders will essentially be higher dimensional. The study of QNMs of these kind of black 
holes projected on the four dimensional brane with the perturbations due to brane localized 
standard model fields were made in [251 ISSl [27] . It was shown that the increase in number 
of extra spatial hidden dimensions dampens the QN frequencies produced via perturbations of 
all kinds of brane localized standard model fields. The presence of charge Q in the black hole 
background also affects the QN spectrum and it is significantly different from the behavior of 
charged black holes in four dimensions. 

However, the study of Dirac QNMs of purely higher dimensional black holes were not 
present in the literature until the work of Cho et al [28]. In that work they have studied 
QNMs of Schwarzschild black holes using the conformal properties of the spinor field. Such an 
idea is perfectly general and in principle can be applied to all higher dimensional spherically 
symmetric black holes. In this paper we will use the abovementioned idea to study QNMs of 
higher dimensional charged black holes arising out of Einstein Hilbert action as well as from 
higher derivative corrections to such actions. In particular we will study the QNMs of higher 
dimensional Reissner- Nordtrom type black holes ^2] and the charged Gauss-Bonnet black 
holes [331 [35l [361 [SZ]- The motivation of the paper is three fold, namely 

To study the Dirac QNM in charged black hole backgrounds in higher dimensions. As 
we have mentioned that study of Dirac perturbations in higher dimension is rather limited, 
this paper will try to fill in a gap in the literature by studying Dirac QNMs in purely higher 
dimensional charged black hole backgrounds in the framework of general relativity and its 
higher derivative corrected scenario. 

To compare the results of charged black hole QNMs in two different scenarios, namely the 
black holes arising out of Einstein Hilbert action and in the higher derivative gravities, precisely 
the Gauss- Bonnet black holes due to Dirac perturbations. We also study the behaviour of the 
QN frequency with the Gauss-Bonnet coupling a for the charged Gauss-Bonnet black hole. 

To compare our results with the available results for brane localized black holes studied by 
Kanti etal [251 [26] and Zhidenko [27]. 
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The plan of the paper is as follows: in the next section we briefly discuss the Reissner 
Nordstrom type solutions and charged Gauss Bonnet solution in dimensions d > 4. In section 
3 we present a brief discussion of WKB method along with a comparative study of the QNMs 
of Reissner-Nordstrom and charged Gauss-Bonnet black holes. Section 4 deals with the quasi- 
normal modes in large multipole number limit and in section 5 we conclude the paper with a 
brief discussion on future directions. Finally in the appendix we give a brief review of the Dirac 
equations in higher dimensional curved background following ^28j. 



2 Charged Black Holes in Higher Dimensions 

In this section we will discuss the charged black holes in higher dimensions. As we have 
mentioned that we will study charged black holes arising out of Einstein's theory of general 
relativity and those arising out of Einstein Gauss Bonnet theory. Let us first discuss the 
Reissner Nordstrom type solutions in higher dimensional Einstein gravity. 



2.1 Reissner-Nordstrom type solutions 

One of the main goals of theoretical physics is to find out a theory which unifies gravity with 
all other fundamental forces in nature. String theory is one of the candidates for such a unified 
theory and it predicts that we live in a world which has dimensions greater than 3 + 1. In this 
context, study of black hole physics is an important area. In this section we discuss higher 
dimensional black holes which are static and spherically symmetric. The study of such black 
holes began with the work of Tangherlini [31] and later by Myers and Perry 



Let us start with the static spherically symmetric metric 

ds' = -f{r)de + r\r)dr' + r'dnl_^, (1) 
The vacuum Einstein equation then implies 

provided that (i > 4. The parameter /i is a constant of integration and is related to the mass 
M of the black hole 

where Ad-2 is the area of the unit {d — 2)-sphere given by 

27r{'^-i)/2 

= (4) 
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One can find the analog of Reissner-Nordstrom solutions in higher dimensions also and there 
the metric is given by [32j 

/=fl-i + ^l (5) 



The electric charge of the black hole is given by 

(rf-2)(rf-3)^2 



^ = -0' (6) 



For 9 < , there is an outermost horizon situated at 



+ ifi' - e'y/' (7) 



This will be of importance for our analysis since we will be interested in the potential just 
outside the outermost horizon of the black hole. In our analysis of the QNM for the charged 
black hole, we will also have to keep in mind the relation ([7j), since this relation gives a constraint 
on the mass and charge of the hole. In other words, one can not work with arbitrary mass and 
charge parameter while working with the Reissner-Nordstrom type metric. The non-extremality 
condition should always be maintained while working with such metrics. 



2.2 Charged Gauss-Bonnet black hole 

We will now briefly discuss the black holes which arise out of Einstein Gauss-Bonnet gravity. 
In space-time dimensions d > 5 the Einstein-Gauss-Bonnet action is given by 

/ = J d^x^R + a' J d'^x^iR^^p^R^"^"' - ARp^R^'< + R^), (8) 

where Gd is the d-dimensional Newton's constant and the parameter a' denotes the Gauss- 
Bonnet coupling. We will choose = 1 from now on and will consider only positive a' which 
is consistent with the string expansion [3i 



The metric for spherically symmetric asymptotically flat Gauss-Bonnet black hole solution 
of mass M is given by Eqn. ([1]), where /(r) has the form 



r2 r2 / 8aM 

where, 

a = lQ'KGd{d-2){d-A)a'. (10) 

For a' > 0, this black hole admits only a single horizon [34] . The horizon r = r/i is determined 
by the real positive solution of the equation 

r^-^ + = 2M. (11) 
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The charged Gauss-Bonnet black hole has the following form of /(r) [37] : 



fir) 



8aM 




AaQ^ 



27r{d - 2){d - 3)r 



2d-4 ' 



(12) 



if M > and a > 0, then there will be a timelike singularity which will be shielded by two 
horizons if Q < Qex- Here Qex is the extremal value of the charge determined from [37j : 



2(d-3) , 

- 3 



2{(i-4) 



Q 



2 

ex 



0, 



where, 



--(d-5)M + 



27r(rf-2)(rf-3) 
^(d 5j 7W +2^(rf_2)(rf-3) 



1/2 



(13) 



(14) 



Again, this equation will be of importance to us since this will constrain the values of the 
mass and charge of the Gauss-Bonnet black hole. For a general discussion on higher derivative 
corrected black holes and their perturbative stability, see Ref [38]. 



3 QNM using WKB method 

Studying QNMs in black hole backgrounds essentially require the solution of Scrodinger like 
wave equation with a particular boundary condition. We will briefly review the derivation of 
the wave equation in higher dimensions due to Dirac perturbations in the appendix. However, 
for most spacetime geometries, the wave equation governing the QNMs is not exactly solvable. 
In our case also the wave equation so obtained is not exactly solvable and we have to look 
for numerical schemes for finding out the QN frequencies. Various numerical schemes have 
been used in the literature to find the QN frequencies, which include direct integration of the 
wave equation in the frequency domain [31], Poschl- Teller approximation [ID], WKB method 
[m 1121 SSI SI] , phase integral method [13 SS] and continued fraction method [U] . We will 
use the WKB method in our case. The WKB method has some advantages over the other 
semianalytic methods since it can be carried systematically to higher orders to improve the 
accuracy. It has been found also that sixth order WKB methods used to find black hole 
QNMs give almost same result as can be found out by full numerical integrations. Also, one 
of the advantages of using the WKB method is that it can give an analytic expression for the 
frequency when one uses the lowest order WKB formula in the large multipole number limit 
which is otherwise almost impossible to find out. 

Having discussed the black hole backgrounds of our interest we now look into the tool that 
we will use in evaluating the QNMs, i.e. the WKB method. As we will see in the appendix, 
the equation we need to solve is 

+ (-^ . V.) f ^ .^F, (15) 
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where is the tortoise coordinate and the potential Vi^2 is given by 



dr^, r 



(16) 



with the choice of outgoing boundary conditions at the horizon and spatial infinity i.e. nothing 
should come in from asymptotic infinity to disturb the system and nothing should come out of 
the horizon. 



(17) 



The formula for QN frequencies using third order WKB approach is given by [H 

= [Vo + (-2Vo")'/'A(n)] - zin + ^){-2V^)'/'[l + Q{n)]. 



where, A = A/i and Cl = Q/{n + ^) and A and Q are given by 



A(n) 



Vr 



(4)' 



8 \ V^' 



(3)' 



(7 + 60z.^ 



288 \ Fo" 



[n 



(3)' 



6912 \ V^' 







384 








' 288 





(77+ 188Z/2) 
1 



(51 + lOOz/^) + 



;i9 + 28z/2) - 



(4)' 



2304 \ V^' 



(67 + 68z/2 



Vr 



(6)' 



288 \ Vo" 



(5 + 4z/2) 



Where, l/J"^ = [d^'V / dx'')^=^^ and z/ = n + 1/2. 

It may be mentioned here that the accuracy of the WKB method depends on the multipole 
number / and the overtone number n. It has been shown [52j that the WKB approach is a good 
one for I > n, i.e. the numerical and the WKB results are in good agreement if / > n, but the 
WKB approach is not so good if / = and not at all applicable for / < n. 



Now, in our case the form of the potential for Reissner-Nordstrom type solutions is given 



by 



dr I r J 



(19) 



K 



A 



r 



{d - l)/xr^ + r" 



A 



where A 



r 



^"^^^ 2/ir'*^^ + 6"^). One can also write A as 

3 3^^^d— 3 ^(i— 3^ 



2d- 



:K-l 



2^6 



{d - 2)e'r 



A 



r 



r 



(20) 
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where r+ and r_ are the outer and inner horizons for the Reissner-Nordstrom black hole in 
ci= dimensions. The explicit form of r± is given by 

T± = {^l±^/J^!^^)^K (21) 

We will need this particular form for the discussion of large angular momentum limit in the 
next section. 

The potential for the Reissner-Nordstrom type black holes and the charged Gauss-Bonnet 
black hole for different dimensions are plotted in figure (1) for a particular value of multipole 
number, charge and Gauss-Bonnet coupling. Since the QN frequencies depend on the height and 
width of the potential, it is clear from the picture, at least qualitatively, that the quasinormal 
frequencies should be different for the two different space times discussed in the paper. Even 
with a very small value of the Gauss Bonnet coupling (a = 0.1 in the figure (1)), there is 
significant changes in the potential. 




(a) Potential for Reissner-Nordstrom black hole (b) Potential for charged Gauss-Bonnet black hole 

with a = 0.1 

Figure 1: The potential for higher dimensional (a) Reissner-Nordstrom and (b) charged Gauss- 
Bonnet black holes. The lower plot is for = 5 and the topmost one is for = 10 for charge 
Q = 0.05 and / = 2 in both cases 

By properly choosing the A one can get the potential for the Reissner-Nordstrom case for 
d = 4, which is given in [20]. Similarly a potential for the charged Gauss Bonnet black hole can 
be found out. Due to the complicated expression for the potential for charged Gauss-Bonnet 
case, we do not explicitly write down the form here. 

Having found out the explicit form for the potential, our next goal is to determine the 
quasinormal frequencies for Reissner-Nordstrom type black holes and charged Gauss-Bonnet 
black holes in higher dimensions. For that we will use the WKB formula for the QN frequencies 
given by Eqn. (|T7|) . In this paper we will be looking at the above mentioned black hole 
backgrounds in space-time dimensions d = 5 to d = 10. 

In figure ([2]) we show the behaviour of the real (2(a)-2(f)) and imaginary part (2(g)-2(l)) of 
the QN frequency u with charge Q for Reissner-Nordstrom type solutions in higher dimensions. 
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Figure 2: The behaviour of Real (a-f) and Imaginary (g-1) part of the quasinormal frequency 
CO with charge Q for Reissner-Nordstrom type black hole in d = 5 to 10 for I = and n = 0. 
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The plots are made for a particular value of the multipole number / = and mode number 
n = 0. As can be seen from the figure the real part of the frequency increases with the increase 
of charge. This implies the real oscillation frequency increases as the charge is increased. 

The behaviour of the imaginary part of the QN frequency uj with charge Q for Reissner- 
Nordstrom type solutions in higher dimensions for / = and n = is shown in Fig. (2(g)-2(l)). 
As can be seen from the figure the negative imaginary part of the frequency falls off with the 
increase of charge. This implies that the damping decreases with the increase of the charge, i.e. 
a longer ringdown phase. For a better and compact look at the behaviour of the frequencies 
with charge and multipole number see Fig. ([3]) 

It may be mentioned here that, in the mathematica code, if we use d = 4, then the values 
for the QN frequencies in four space time dimensions can be found out and our result matches 
with the result obtained by [20] and [21]. Both the authors of [Sni El] have used Poschl Teller 
approximation scheme to determine the QN frequencies whereas we have used WKB method 
to determine the same. We have also checked our results with the results obtained by Cho et 
al in [28] by explicitly putting Q = in the metric and have found that the results matches 
with the results of higher dimensional Schwarzschild black holes as obtained by them. 



Figure 3: Behaviour of the (a) Real and (b) Imaginary part of u with / and Q for the higher 
dimensional RN black hole 

Let us now discuss the behaviour of the real and imaginary part of the QN frequency 
with the space time dimensions. The real part of the frequency increases with the increase of 
dimension while the imaginary part also increases with the space time dimension (See Fig. H]). 
So, the real oscillation frequencies are increasing as the space time dimensions increase, but the 
damping also becomes larger for higher dimensional charged black holes. Thus it is found that 
the qualitative behaviour of the QN frequencies for purely higher dimensional RN black hole 
is the same with the behaviour of QN frequencies for higher dimensional Reissner-Nordstrom 
black holes projected on the brane [25]. Though it may be mentioned that the quantitative 





(a) Variation of Re uj with Q and I 



(b) Variation of Im lu with Q and I 
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Space time Dimension 



Space time dimension 



(a) Variation of Re lo with d 



(b) Variation of Im lo with d 



Figure 4: Behaviour of the (a) Real and (b) Imaginary part of uj with space time dimension for 
/ = 2 and n = for the higher dimensional RN black hole 
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Figure 5: Re w vs Im a; for different I for the higher dimensional RN black hole. Each four 
segments in the plots correspond to / = 0, 1, 2, 3 respectively. 
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nature of the frequencies are hugely different. Both the real and imaginary parts of the pure 
d-dimensional RN black holes are larger in magnitude when compared with the results of [25] , 
where the QN frequencies for RN black holes projected on the brane were calculated. We have 
checked this for dimensions d = 5,6 (for which the data were available in [23]). 

The real and imaginary part of the quasinormal frequency is plotted in figure (jS]). Each four 
distinct segments in all the plots correspond to multipole values 1 = 0, 1,2 and 3 respectively. 
The points in each segment correspond to different values of charge starting from Q = to 
Qex, where Q^x is the extremal value of the charge. 

Now, let us take a look at the results for the charged Gauss-Bonnet black hole. Due to 
a huge set of data, we do not give the data for all the dimensions with all values of charge, 
coupling and multipole number in this paper, however we only consider a particular d = 7 and 
a particular multipole number / = 1 and overtone number n = for the sake of simplicity. 
The study of scalar field evolution in the Gauss Bonnet background was done in [H] and [19] 
while vector and tensor perturbation was done in [50] (for recent study see [5T]). It may be 
noted that there were no such studies on QNMs in the charged Gauss-Bonnet background due 
to Dirac perturbations, this paper tries to fill in this gap in the literature. 

As can be seen from the Figs. ([ni^a)-[SI^f)), the real part of the frequency for the charged 
Gauss-Bonnet black hole increases with the increase of charge whereas Figs. ([6](g)-l6](l)) suggests 
that the imaginary part decreases with the charge implying that the real oscillation frequency 
will be increasing with charge whereas the damping decreases with it. It is in this sense that 
the behavior of the frequencies in both the Reissner-Nordstrom and charged Gauss-Bonnet 
background is same. Now let us comment on the Gauss-Bonnet coupling a — »• limit. It has 
been shown in [48] and [50] that in the a ^ limit, the quasinormal modes for scalar field 
perturbations and vector and tensorial perturbations of the uncharged Gauss Bonnet black 
hole yields Schwarzschild QN frequencies. This is very easy to understand since in the limit 
a ^ 0, the Gauss Bonnet metric looks like (1 - 2M/r'^-^ + 4aMVr^'^"^ + ■■■)• The third 
term in the above metric is actually O(a^) term and hence for very small values of a the 
quasinormal frequencies for Gauss-Bonnet black hole goes to the Schwarzschild values. The 
same thing is expected here because in the a ^ limit, the charged Gauss-Bonnet metric 
looks like (1 - 2M/r'^-^ + Q^/r^-^-^ + 0{a^) + ■■■), where Q = Q / ^2TT{d - 2){d - 3). We have 
checked for different small values of a (a = 0.0001, 0.0002, ■ ■ ■ , 0.0005), keeping the charge fixed 
that the quasinormal frequencies for Reissner Nordstrom black holes in different dimensions are 
indeed produced. 

Having discussed the behaviour of the real and imaginary parts of the frequencies with 
the charge of the Gauss Bonnet black hole, let us now concentrate on the variation of the 
frequencies with the parameter a, which is the Gauss-Bonnet coupling. The variation of the 
real and imaginary parts of the frequency with the Gauss-Bonnet coupling is plotted in Fig. 
(I7j). As can be seen from the figures, the real part of the frequency increases with the increase of 
Gauss-Bonnet coupling while the negative imaginary part decreases with a. It may be noticed 
that the variation of the imaginary part in higher dimensions is comparatively slower than 
the lower dimensions. For example, in c? = 8, 9, 10, the negative imaginary part first starts 



12 



0.00 0.02 0.04 0.06 0.08 



0.00 0.02 0.04 0.06 O.OS 




(a) Re{uj) vs Q /or a = 0.1 (b) i?e(w) vs Q for 0.2 (c) Re{uj) vs Q for 0.3 



0.00 0.02 0.04 0.06 0.08 



0.00 0.02 0.04 0.06 O.OS 




0.00 0.02 0.04 0.06 O.OS 



(d) i?e(w) vs Q for a = 0.4 (e) Re{uj) vs Q for a = 0.5 



0.7388 - 

0.7387 - 

0.7386 : 

0.7385 : 

0.7384 : 

0.7383 : 

0,7382 : 



0.00 0.02 0.04 0.06 



0.00 0.02 



0.04 0.06 
Q 



(f) Re{uj) vs Q for a =1.0 



0.02 0.04 0,06 O.OS 

Q 



(g) -'Im{ijj) vs Q for a = 0.1 (h) —Im{uj) vs Q for a = 0.2 (i) —Im[ijj) vs Q for a = 0.3 



0.6798 . 
0.6796 - 



0.00 0.02 



0.04 0.06 



0.00 0.02 



0.04 0,06 
Q 



0,02 0,04 0,06 

Q 



(j) —Im{uj) vs Q for a = 0.4 (k) —Im{uj) vs Q for a = 0.5 (1) —Im[ijj) vs Q for a = 1.0 



Figure 6: The behaviour of Real (a-f) and Imaginary (g-1) part of the quasinormal frequency 
uj with charge Q for charged Gauss-Bonnet black hole in = 7 for / = 1 and n = 0. 
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decreasing with a, then remains unchanged for certain consecutive values and again starts 
decreasing with a. 
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(a) Plot of Real w vs a for charged (b) Plot of Imaginary w vs a for (c) Plot of Quality Factor vs a in 
Gauss-Bonnet black hole charged Gauss-Bonnet black hole d = l 

Figure 7: The variation of (a) real, (b) imaginary parts of frequency and (c) quality factor with 
Gauss-Bonnet coupling a for a fixed value of charge Q = 0.05 and multipole number / = 1. 
Dimension of the space time increases for (a) and (b) from bottom to top, i.e. the lowest plot 
in each graph corresponds to ci = 6 and top one corresponds to = 10 

If one defines the quality factor as Q ~ | ^^^^ | , then its variation is plotted in fig ([7t). 
This shows that in spite of the behaviour of the real and imaginary parts of the frequency 
stated above, the quality factor increases as one increase the value of a. It is also interesting 
to note that though Q increases faster for small values of a, it becomes somewhat saturated 
as a is gradually increased, which is similar to the behaviour of the quality factor observed for 
uncharged Gauss-Bonnet black hole on the brane also [27] . 



4 QNM in large multipole number limit 

Let us now focus on the large angular momentum (multipole number) limit, i.e. now n — > oo. 
Our aim in this section is to simplify the otherwise complicated expression for the frequency 
given by Eqn. (fT7|) and get an analytic expression for u. For this, we will consider the large 
multipole number limit, where one can see from the expression of the potential ( fTOll . the most 
dominant term of the potential in k — > cxd limit will be 

V^iU^oo-/^. (22) 

Hence this approximation simplifies the potential a lot. Next, we will consider the WKB formula 
only upto the first order, i.e. now we will write Eqn. (fTTll as 

u;'^Vo-^(n + ^^ {-2V^)'/' + ■■■, (23) 
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where Vq is the maxima of the potential Vi that we are working with. For finding the maxima 

dVi 
dr 



of the potential we solve the equation ^ = and get 



1 

2d — 6\ "^-^ r>r_ 



(24) 



d-3 _ 



The form of the potential with this maxima of r is 

V^il™ = (^2^r_r+ (r^-3 + r^-3) 3^j 4-2d 

((2^r_r+ + ^^-3) 3^ j d-3 _ ^d-3j ^35) 

Next one uses this to find out an expression for the frequency 

= [2'^r_r+{r'^~^ + r'f_"^) 3^]^"^'^ x [(2^r_r+(r'^"^ + r^|_"^) 3^)''"^ - r'^"^] x 
[(2^r_r+(r^-3 + r^-3)3^)^-3 _ ^^-3] 

— z n + - X 
4 V 2; 

^ ^ ^_2l^(rl-3 + r:f.-3)-dr3(2d^r„r+(rl"^ + r^"^)^)"2'^(-(c/ - l)d(2^r_r+ 



^2 ^2 



(2^r_r+(/-3 + r^-3)3^)2rf(^d-3 + ^^-3^^ ^ 2^{d{2d - 7) + 6)r^+Vf =^)]5 (26) 

Which indeed is very complicated, though we have simplified our potential and used the 
first order WKB formula. We can note from the above form that the QN frequencies indeed 
depend on the mass and charge of the black hole via r+ and r_ (since r+ and r_ are determined 
in terms of /i and 6) . One can check by putting d = 4 in the above formula that the result for 
4 dimensional Reissner Nordstrom black hole is indeed obtained. 



5 Conclusion 

In this paper, we have discussed the Dirac quasinormal frequencies of charged black holes 
in higher dimensions. We have investigated two different backgrounds, namely the higher 
dimensional Reissner-Nordstrom and the charged Gauss-Bonnet black holes. Though the field 
of studying QNMs is being saturated day by day, a comparative study of QNMs in different 
scenarios might be interesting. 

One of the main ideas of [28J was to study fermion quasinormal modes in purely higher 
dimensional Schwarzschild background within the framework of split fermion models where the 
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quarks and the leptons are forced to live on separate branes in order to keep proton stability. 
The present work can also be interpreted as studying split fermion quasinormal modes for 
higher dimensional charged black holes. Within this model, we found that the real part of the 
QN frequencies increases with the increase of the charge in both the backgrounds discussed in 
this paper, while the imaginary part decreases with the charge for both. The behavior of the 
frequencies are studied with the increases in space time dimensions. We found that the real 
part of the frequency increases as we increase the space time dimensions while the damping 
also increases. A comparative study of the quasinormal modes in the higher dimensional RN 
and charged Gauss-Bonnet black hole suggests that the variation of the frequencies with charge 
is much more rapid in the Reissner-Nordstrom background than the charged Gauss-Bonnet 
background. Quantitatively we found that the real and imaginary parts of the frequencies 
in the split fermion models are larger than their brane localized partners. One can also get 
back the QN frequencies of the black holes arising out of pure Einstein theory of gravity from 
the values of the QN frequencies of charged Gauss-Bonnet background in the Gauss-Bonnet 
coupling a — i> limit. We have also studied the behaviour of the real and imaginary parts of 
the frequencies with the Gauss-Bonnet coupling a and found that the real part increases with 
a, while the negative imaginary part decreases. For higher dimensions this variation in the 
imaginary part is slower and it can be seen that for some consecutive values of a, the change 
in the imaginary part is very small and then if we vary a again, the negative imaginary part 
starts decreasing again. 

It would be interesting now to calculate the black hole absorption cross section for bulk RN 
and Gauss Bonnet fermions in higher dimensions and a comparative study of late time fall off 
would also be interesting in these backgrounds. As it was found that the massive fields might 
change the low lying modes of the QN spectrum, the study of massive Dirac perturbations in 
these background will also be an important aspect. 
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7 Appendix A: Dirac Equation in spherically symmetric 
higher dimensional black hole background 

In this section we will study the Dirac equation in higher dimensions in a static spherically 
symmetric black hole background. Essentially this section contains a brief review of the works 
done in [281. 
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Let us start with a d dimensional spherically symmetric metric of the form 



2 

d-2 5 



ds'^ = -f{r)dt^ + f-\r)dr^ + r^dn 

where dQ'^_2 is the metric for {d — 2) sphere. Following 
transformation under which the metric and the determinant of the metric behaves as 



(27) 



301. let us think of a conformal 



9 



det Ifif^^l 



^'^9^^,y 



-9 



(28) 



where Q is the conformal factor. Now, under such a conformal transformation the spinor field 
behaves as ip ^ ip = where n can be determined by claiming that the Dirac Lagrangian 
remains invariant under such a conformal transformation. Using the fact that {7^,7'^} = 2g'^'' 
and 7^ = f2~^7^ along with the massless Dirac Lagrangian one can get n = —{d — l)/2. Thus 



Following [28], we choose 



then = 1 



(i-d) 

(d+l) 



(29) 



'^^/'^ip. The conformal metric ds'^ then has a 



completely separated t — r part and the {d — 2)-sphere part. The massless Dirac equation then 
can be written as 



r^n^ = 



(30) 



where V, 



lVac^lt,^''\ and 



~ab 



* [7°, 7^]. The problem now 

is to find out the 7 matrices for higher dimensional space time. In [53], a complete discussion 
about this problem was given for flat space time. Following their notation and using [28] we 
write the separated Dirac equation as 



-a 



(7*V 



iJSd 



ij = o 



(31) 



Where, cr* are the Pauli matrices 



a 



1 

1 



a 



-i 

1 



a 



1 

-1 



(32) 



Following 



one can write 



±i I 



d-2 



0, 1, 2,-- 



(33) 
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where xf the eigenspinors for the {d — 2)-dimensional sphere. Again, following |28], one 

can write = ''^^{4>ti''^y^)xt + 0r(^5'^)xr) using the orthogonality of the eigenspinors. The 

I 

Dirac Equation can then be written as 

0+ - ^CT' (/ + ^) 0+ = za' (34) 

In writing the above equation we have changed the notation by removing the tilde from the 
expressions and also we are working here with the positive sign of the eigenspinors. It may be 
mentioned that one can also choose to work with the negative sign as well. Next, we use the 
ansatz 

The Dirac equation then simplifies to 



V7 



d Nf 



2^dr 



+ (36) 



dr r \ 2 
+ + (37) 



dr r \ 2 

By defining the tortoise coordinate as dr^ = dr/f{r) and introducing a function 

W = —K (38) 

r 

where k = (/ + (c/ — 2)/2), one can finally arrive at 

'^^ + V^G = u:'G- (-^ + V2)f = u'F, (39) 



dr"^ J ' \ dr^ 

where 

V,,2 = ±^ + W\ (40) 

The potentials Vi and V2 corresponding to Dirac particles and anti-particles are supersymmetric 
to each other and dervide from the same superpotential W. It has also been proved |21j that 
the Dirac particles and anti-particles have the same QN spectra. In this paper we have used 
only Vi to study the quasinormal modes. 
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